Phase-space for the breakdown of the quantum Hall effect in epitaxial graphene 
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We report the phase-space defined by the quantum Hall effect breakdown in polymer gated epi- 
taxial graphene on SiC (SiC/G) as a function of temperature, current, carrier density, and magnetic 
fields up to 30T. At 2K breakdown currents (Ic) almost two orders of magnitude greater than in 
GaAs devices are observed. The phase boundary of the dissipationless state {pxx = 0) shows a 
{1-{T/Tc)^) dependence and persists up to Tc > 45K at 29T. With magnetic field Ic was found 
to increase oc B'^'^ and Tc oc B^'^^. As the Fermi energy approaches the Dirac point, the u — 2 
quantized Hall plateau appears continuously from fields as low as IT up to at least 19T due to a 
strong magnetic field dependence of the carrier density. 
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The quantum Hall effect (QHE) observed in two- 
dimensional electron gases is defined by vanishing lon- 
gitudinal resistivity, p^x = 0, and a quantised Hall resis- 
tance, pxy = h/ve^ for v =integer. Ever since its first 
observation [1] in silicon the QHE has been used as a 
quantum electrical resistance standard which has been 
most extensively developed using GaAs devices [5]. In re- 
cent years since the the first isolation of graphene and 
the observation of the integer QHE "3] [3] the attention 
of quantum Hall metrology labs has turned to graphene 
as potentially a more readily accessible resistance stan- 
dard capable of operating at higher temperatures and 
measurement currents with lower magnetic fields. This 
is in part due to its large cyclotron energy gaps aris- 
ing from the high electron velocity at the Dirac point. 
Recent experimental workj^ has also shown that it has 
high electron-phonon energy relaxation rates, an order of 
magnitude faster than in GaAs heterostructures, which 
play an important role in determining the high current 
breakdown of the QHE. In particular polymer gated epi- 
taxial graphene on SiC has been shown to be an excep- 
tional candidate for metrology lU E] and the universality 
of quantisation between it and GaAs has been shown to 
be accurate within a relative uncertainty of 8.6 x 10"^"'^ [S]. 

If epitaxial graphene is to be used as a quantum re- 
sistance standard it is important to understand the ex- 
perimental limits which confine the phase-space where 
the accurate, dissipationless QHE can be observed. Such 
a phase-space is determined by temperature, T, carrier 
density, n, magnetic field, B, and current /. The break- 
down of the QHE is defined as the point where devia- 
tions from quantisation, lS.pxy can be observed and this 
is strongly correlated with the point where pxx 7^ 0. A 
linear relationship of lS.pxy/ pxy oc pxx is typically ob- 
served in GaAsfS] and recently in graphene[IOj. there- 



fore measurement of the / — Vxx characteristics in the 
quantum Hall regime also determines the maximum cur- 
rent consistent with maintaining a quantised pxy. At 
high currents a sudden onset of longitudinal resistance 
is observed above a critical current /c[2l 13 HI]- Mod- 
elling the breakdown using the bootstrap electron heat- 
ing (BSEH) model from Komiyama and Kawaguchi [T^ 
which looks at the balance between the rate of increase 
of energy input and electron energy loss has been rea- 
sonably successful. The Quantum Hall effect has already 
been reported in graphene at room temperature using 
magnetic fields of 45T[TT, however the plateaus did not 
show exact quantization as the resistivity was still finite 
(^ \QVt). In this work we address the formation of the 
zero-resistance state which corresponds to the dissipa- 
tionless current flow needed to maintain the full quantum 
Hall condition. 

Two devices were studied, prepared from epitaxially 
grown graphene on the Si-terminated face of SiC. Each 
device was lithographed using an e-beam and oxygen 
plasma etching into an eight leg Hall bar geometry 
{W/L = 4.5) with widths of I^ = bpui and W = 35^m 
for Sample 1 and Sample 2 respectively. Samples were 
electrically connected with large area Ti-Au contacting. 
A polymer gating technique using room temperature UV 
illumination was used to vary the electron density from 
1 — 16 X lO^^cm^^ as described in Ref.[T4]. D.C. mag- 
netotransport and I — V data was taken using magnetic 
fields from a 21T superconducting solenoid and a SOT 
20MW resistive-coil magnet at the LNCMI Grenoble. 

Fig. [Tb shows Pxx and pxy for Sample 1 with nB=o — 
6.5 X lOT'^cm"^. We observe Shubnikov-de Haas oscilla- 
tions in filling factors up to z/ = 8, and a j/ = 2 quantum 
Hall plateau beginning at B — 8T with pxx = from 
B = lOT. This 1/ = 2 state is over 20T wide and ob- 
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FIG. 1: a) / — Vxx characteristics of Sample 1 at 2. OK, with 
a breakdown condition of Vxx ~ 10^J.V, giving a maximum 
critical current density j^ ~ 43A/m at 23T. b) Combined 
magneto-transport {pxy{h\ue) ,pxx{hla,ck)) data and /— Vxx — 
B contour plot, the hashed region represents Vxx < lO^V, 
the dissipationless quantum Hall regime, c) Magnetic field 
dependence of the carrier density (black), following lines of 
constant filling factor (red) while Ep lies between Landau 
levels and then the charge transferred from surface donors in 
SiC, n{B,N) (green), while the Landau levels fill, from the 
model in [?]■ 



servable all the way up to the maximum magnetic field 
of 30T. A series of / — Vxx traces were taken every Tesla 
along the plateau to investigate the breakdown, with typ- 
ical examples in Fig. [l}i at T = 2K. At 23T we find 
Vxx — until I = Ic = 215/iA, where we define the criti- 
cal breakdown current at Vxx(Ic) = 10/iV, just above the 
noise-level of our measurements (Fig. Wp.), correspond- 
ing to a resistivity of pxx ~ 0.01J7. Such a high break- 
down current for a device just 5/ini wide, giving a crit- 
ical current density of jc = AZA/m, is truly exceptional 
in comparison to even the most well optimised GaAs de- 
vices (jf "-^^ -l-2A/m). The full set of / - Vxx traces 
are plotted in Fig. [TJd as a contour plot. The hashed 
region is the the phase space where the dissipationless 
QHE is observed. The critical current Ic increases along 
the plateau with a peak around 23T. Unlike traditional 
semiconductor quantum Hall systems which show a very 
sharp peak in Ic centered at integer filling factor [IS], 
the peak breakdown current occurs at fields much greater 
than V = 2 calculated from the zero-field carrier density 
and changes very little in magnitude over a wide range 
of fields. This is due to the strong magnetic field depen- 
dence of the carrier density in epitaxial graphene grown 



on Si-terminated SiC [TB]. Carriers are transferred to 
the graphene from the surface donor states of the SiC 
which are assumed to have a constant density of states. 
The charge transfer, ns{B,N) is proportional to the dif- 
ference between the workfunction of the graphene and 
the SiC. This causes the unbroadened Landau levels to 
be completely filled over a wide range of magnetic fields 
[7], particularly when the Fermi Energy Ep is between 
the N=0 and N=l Landau levels as in the region above 
IIT in Fig. [T]:. Assuming that the peak Ic occurs at 
V — 2 suggests that the carrier density has increased to 
n = 1.1 X lO^^cm"^ by 23T and is still increasing. As a 
result the breakdown current is relatively independent of 
magnetic field which adds to the convenience of epitaxial 
graphene as an electrical resistance standard. 

At the lowest carrier density studied using Sample 2 
(nB=o -- 1 X 10"cm-2), \he v ^ 2 state (Fig. [2J1) be- 
gins at B =1T and persists up to the maximum field 
studied for this sample of 19T. The breakdown current 
shown in Fig. [2k is negligible at low fields (B <3T) but 
rapidly increases reaching a peak at B =7T, suggest- 
ing a carrier density of nB=7T = 3.5 x 10^^cm~^. At 
7T Ic = 140/iA, giving jc = 4A/m for this 35/i-m wide 
device. Importantly, from an applications perspective, 
Ic ~ 100/iA by 5T, a magnetic field which is readily ac- 
cessible with simple benchtop magnets. Applying the 
charge transfer model [TB] the magnetic field for peak 
breakdown is accurately predicted (fig. |2|d) but above 
this no further increase in carrier density is expected due 
to the finite density of donor states. The data suggest 
that the carrier density is still increasing, as the break- 
down current has only decreased by a factor of 0.55 by 
19T, probably due to the influence of level broadening 
which is not included in the original model [TH]. In typ- 
ical semiconductor 2DEGs [H [15], breakdown currents 
show a triangular behaviour with a plateau width (de- 
fined by Ic{v)/Ic ^ 0) of Av/v ^ ±0.2. Assuming a 
level degeneracy (77) of 4 for the v — 2 plateau due to 
the valley and spin degeneracies in graphene, this should 
correspond to a total plateau width of Az/ = ±0.8, and 
Ic should halve by 9T (y = 1.6). This is consistent with 
results reported for exfoliated graphene [T7] . By contrast 
the slow decrease in Ic seen in Figj2^ suggests that the 
occupancy remains u ^ 1.6 up to 19T where the carrier 
density has increased to n ^ 7 x lO^^cm^^. 

An Arrhenius analysis of the activated conductivity at 
higher temperatures (50-80K), above the variable range 
hopping regime[18l [19], was used to estimate the mag- 
netic field dependence of the Fermi energy, Ep, by mea- 
suring the activation gap A as a function of magnetic 
field. We assume that this measures the separation of Ep 
from the conducting states E^ of the nearest Landau level 
(N=l for B< 7T, N=0 for B> 7T), where A = \Ef,-Ep\. 
Fig. [2j;, shows A and the value of Ep which has been 
deduced by assuming that it is midway between the two 
Landau levels at 7T where v = 2. At low fields Ep 



corresponds to the approximately constant value of 40 
meV deduced from the low field carrier density. Above 
2.5T the carrier density begins to increase due to charge 
transfer from the substrate which keeps the Fermi energy 
in the gap between N=l and N=0, and the system en- 
ters the dissipationless quantum Hall state. Above 7T 
Ep falls slightly but appears pinned close to a constant 
energy of Ep ~ 40 meV suggesting that there may be 
a specific surface impurity level close to this value. This 
suggests that the i/ = 2 plateau could extend up to higher 
fields still until the extended states of the symmetry bro- 
ken N=0 state pass through the pinned Fermi level. 

In addition to the high currents observed above, high 
temperature (T > AK) operation is required for an acces- 
sible resistance standard. The temperature dependence 
of the breakdown was studied at several carrier densities 
for the peak Ic a.t i/ = 2, and for the highest carrier den- 
sity of n = 1.6 X lO^^cm^^, at 29T as the maximum Ic 
was just beyond our maximum field. In GaAs and InSb 
the temperature dependence [TSl HOI IS] of the critical 
current has been shown to be of the form 



/e(r) = /,(0) 1 - 
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where Tc is the temperature at which Ic — 0. Fig. 
[3^ shows that this equation also describes the temper- 
ature dependence of Ic very well in epitaxial graphene 
for the magnetic fields studied. Rigal et al. [H] justified 
this based upon phenomenological similarities with the 
Gorter-Casimir two-fluid model for superconductors, and 
Tanaka et al. |20] have proposed a model which predicts 
this behaviour based on a temperature-dependent mobil- 
ity edge caused by the temperature dependence of the 
tunneling probabilites from localised to extended states 
at the centre of the Landau levels. Experimentally only 
limited evidence exists for the dependence of Tc on mag- 
netic field with values for GaAs[20l [21] in the range 4.8 
- 7.7K for B values of 4.8 - 7.7T at i/ = 4 and Tc ~ 1/z/, 
while for InSb Tc = 8K at 6.1T[T3], suggesting a linear de- 
pendence of Tc upon field. It is therefore surprising that 
for graphene we see a strong superlinear scaling, as shown 
in Fig. |3p with a best fit of Tc oc B^-^^, which extrapo- 
lates to Tc = 105K at 45T. The rate of increase of the the 
cyclotron energy gap between the A^=l and A^=0 Landau 
level is sublinear, given by En = sgn{N) x c* ^/2ehB\N\, 
where c* is the electron velocity, suggesting a weaker 
overall field dependence. One significant difference in epi- 
taxial graphene is the magnitude of the disorder which 
means that the activation energy A — \E^ — Ep\ at 1^—2 
has a large offset due to level broadening and is known to 
increase more rapidly than the cyclotron energy [TH] due 
possibly to smaller broadening for the N=0 Landau level 
which is topologically protected|3j. 

By contrast Ic has been extensively studied and is well 
known experimentally to scale as B^/^[31[TT1 [50] as pre- 
dicted by several of the models for breakdown[T2] which 
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FIG. 2: a) Magnetotransport {pxy{hlue) and p^^a; (black)) and 
corresponding / — Vxx — B contour plot for the 35/ini wide 
device at T =L5K. b) Theoretical prediction of magnetic field 
dependent carrier density after Ref. !16i. c) A as a function 
of magnetic field and the resulting Ep ■ 
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TABLE I: Material comparison for QHE breakdown at i^ = 2 



include factors for the cyclotron energy and the inverse 
magnetic length. Fig. [st shows the values for jc—Ic/^ at 
I' = 2 for both samples after each UV illumination. The 
highest values observed are also consistent with a B^^^ 
dependence, although there is considerable spread in the 
experimental values, probably because the extended UV 
illumination is thought to introduce significant spatial in- 
homogencities which are likely to reduce jc- Interestingly 
despite the spread of jc values the same samples produced 
the very clear systematic dependence of Tc shown in Fig. 
[3J It should be noted that the values are significantly 
higher for the 5/Lim Hall bars and there is some evidence 
that quantum Hall breakdown current densities are larger 
for smaller Hall bar widths [^[^ 

The most widely accepted theory for the QHE break- 
down is the bootstrap electron heating model proposed 
by Komiyama and Kawaguchi|12| in which the quantum 
Hall state becomes thermally unstable above a critical 
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FIG. 3: a) Normalised temperature dependence of breakdown 
current at several magnetic fields, fitted with equation [T] (b) 
The magnetic field dependence of Tc, with a best fit of 0.08 * 
B^'*^. (c) The magnetic field dependence of jc for the two 
Hall bars with the sequence of illumination shown by arrows. 



Hall electric field where the rate of change of electron- 
phonon energy loss rate becomes less than the rate of 
increase of input power. This predicts a critical break- 
down electric field of 
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where t^. is a characteristic electron-phonon energy re- 
laxation time. Recently much experimental [531 [55H27] 
and theoretical [28j interest has focused on the way hot 
electrons lose energy to the lattice in graphene. We use 
the values of Tg, observed at Tc from the damping of 
Shubnikov-deHaas oscillations [HI |131 HHIj to calculate 
the predicted Ic ioi u — 2 and compare these to conven- 
tional semiconductor 2DEGs in Table [ij The graphene 
values are considerably larger, as compared for example 
to InSb, which has the lowest mass of the III-V semicon- 
ductors, m* = 0.02r?T,e[5D]. At 7T the cyclotron energy 
gap is 105meV for graphene, compared to 40nieV in InSb, 
however we find an order of magnitude increase in current 
density for graphene over InSb. This is mainly a result 
of the factor 6 difference in Tg between the two systems. 
The increase of Tc with field causes t^ to decrease and 
the dependence of jc on magnetic field to be superlinear. 
In summary, we have investigated the phase space in 
which the dissipationless quantum Hall state exists for 
epitaxial graphene. The data support the idea that this 
system can be described in terms of a phase diagram 
where the temperature dependence of the critical cur- 
rent follows a ex (1 — (^)^) behaviour as seen in GaAs 
and InSb quantum Hall systems. We demonstrate that 
both the critical temperature and current are strongly 
magnetic field dependent and that at high fields criti- 
cal current densities can be more than a factor 30 larger 



than previously observed in other systems. In epitax- 
ial graphene, charge transfer from the carbon layer be- 
tween the graphene and the SiC substrate also leads to a 
strongly magnetic field dependent carrier density and an 
exceptionally wide i/ — 2 plateau due to charge transfer 
from surface impurities followed by pinning to a constant 
energy associated with a surface impurity level. 
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